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Abstract. By exploiting density-matrix renormalization group techniques, we investigate the 
spin dynamics of a spin-1/2 one-dimensional J 1 -J 2 XXZ model with competing ferromagnetic Ji 
and antiferromagnetic J 2 exchange couplings under applied magnetic fields. Numerical results of 
spin excitation spectra show that in the field-induced spin quadrupole regime, the longitudinal 
component has a gapless mode and the transverse component has a gapped mode irrespective 
of the exchange anisotropy. The excitation gap of the transverse spin excitation increases as 
the exchange anisotropy increases over the XY-like and Ising-like regions, demonstrating that 
two-magnon bound states are stabilized due to the easy-axis anisotropy. 


1. Introduction 

Frustrated quantum spin systems in low dimensions have provided a fascinating playground to 
explore strange new phases that have no classical analogue, and rich physics of quantum phase 
transitions between them. One of the simplest models with frustration and low-dimensionality 
is a spin-1/2 one-dimensional J1-J2 Heisenberg model, in which the nearest-neighbor exchange 
interaction J\ and the next-nearest-neighbor exchange interaction J 2 compete with each other. 
Note that when J 2 is antiferromagnetic, frustration occurs in either case of antiferromagnetic or 
ferromagnetic J\. Recent studies on the ferromagnetic J\ case in a magnetic field have revealed 
that the ground state is a vector chiral (VC) state in low magnetic fields, while it turns to an 
exotic Tomonaga-Luttinger-liquid state in high magnetic fields mmmm- In the high-field phase, 
the formation of two-magnon bound states leads to a spin quadrupole state. The longitudinal 
spin and quadrupole correlations are quasi-long-ranged, while the transverse spin correlation is 
short-ranged. The so-called spin nematic (SN) state is realized in the vicinity of the saturation, 
while the ground state shows a crossover to a spin-density-wave (SDW) state with lowering the 
magnetic field. The ground-state phase diagram in the presence of the exchange anisotropy has 
also been analyzed numerically |5j. As for the experimental realization of the quadrupole state, 
LiCuV 04 has been studied as a possible candidate material. Neutron diffraction experiments 
have shown that the field dependence of an incommensurate magnetic peak appearing in high 
magnetic fields is well explained by theoretical predictions of the quadrupole state [6]. Inelastic 
neutron scattering measurements have been carried out at zero magnetic field, indicating the 
existence of a two-spinon continuum [?], whereas the excitation spectrum in the high-field phase 
has not been reported yet. 

To clarify the property of the field-induced quadrupole state from the viewpoint of the spin 


dynamics, we have studied the dynamical spin structure factor of the spin-1/2 one-dimensional 
J\-J -2 Heisenberg model by exploiting numerical methods such as a dynamical density-matrix 
renormalization group [8]. In this paper, to gain further insight into the effect of the exchange 
anisotropy, we analyze the J\-Jo XXZ model. We present numerical results on the anisotropic 
behavior of longitudinal and transverse spin excitation spectra, and discuss the dependence on 
the exchange anisotropy. 


2. Model and method 

We consider a spin-1/2 Ji-J 2 XXZ model on a one-dimensional chain with N sites, described by 


H = Ji £(S?Sf +1 + s?sv +l + A SfSf +1 ) + J 2 £( 5 ? 5 ?+ 2 + ^53+2 + ^SfSf +2 ) - h ]T Sf, ( 1 ) 


where Sf is the a component of the spin-1/2 operator at site i, J\ (< 0) and J 2 (> 0) denote 
the ferromagnetic nearest-neighbor and the antiferromagnetic next-nearest-neighbor exchange 
interactions, respectively, A is the anisotropy of the exchange interaction, and h is the magnetic 
field along the z direction. In the present study, we fix J\ = — 1 and J 2 = 1, and examine the 
dependence on A. 

We investigate dynamical properties by using density-matrix renormalization group (DMRG) 
methods in open boundary conditions mm- To clarify the characteristics of anisotropic 
magnetic excitations, we calculate the z and x components of the dynamical spin structure 
factor, given by 


S a (q,u) = - lim-Im^clffi , F 1 ,, , ■ Sq IVm), (2) 

7T W + E G — H + 1?7 H 

where \i/jg) is the ground state and Eq is its eigenenergy. We set an infinitesimal value q to 0.1. 
The momentum representation of the spin operator in open boundary conditions is given by 


s q = 55 sin (^ : ) • ( 3 ) 

For the calculation of the dynamical spin structure factor at a given magnetization m = M/N 
with M = Sf, we set the magnetic field h to be the midpoint of the magnetization plateau 
of m. We employ a dynamical DMRG method m, in which we need three target states, i.e., 
|-0g) J Sq IV’g)) an d [cj + Eq — H + irj\ _1 5'"I^g)- We perform dynamical DMRG calculations with 
32 sites at m = 0 and m = 0.25 for typical values of A. 

3. Numerical results 

First, let us look at spin excitation spectra at zero magnetic held. In Fig. UJ we present intensity 
plots of S z (q,uo) and S x (q,u ) at rn = 0. Figured! a) shows S z (q, ui) at A = 1 for the isotropic 
Heisenberg case, where the system is in a Haldane dimer phase m- Here, S z (q,co) = S x (q,ui ) 
due to the SU(2) rotational invariance, as shown in Figs. DJa) anddlb). There is a sharp peak 
at (qo,uo) = ( 167 t / 33 , 0.14), and a large amount of spectral weight is concentrated on this peak. 
We note that coq seems to be a finite value because of the finite-size effect. In fact, wo becomes 
small as the system size increases, and it is extrapolated to almost zero in the thermodynamic 
limit. This is consistent with the renormalization group study suggesting an exponentially small 
gap [12]. The lower boundary of a continuum is found to form a sinusoidal dispersion, which 
corresponds to the des Cloizeaux-Pearson mode of spinon excitations in the limit of J 1 /J 2 = 0, 
i.e., decoupled antiferromagnetic chains m- The sinusoidal dispersion is clearly visible due to 
the large intensity for q < qo, while that for q > qo is less distinct, since the spectral weight is 
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Figure 1. The dynamical spin structure factors S z (q, u) and S x (q,u ) at J\ = —1, J 2 = 1, and 
m = 0 for typical values of A: (a),(b) A = 1 for the isotropic Heisenberg case; (c), (d) A = 0.6 
for the XY-like case; and (e), (f) A = 1.4 for the Ising-like case. Three panels in the left column 
[(a), (c), (e)] are S z (q,oj), and those in the right column [(b), (d), (f)] are S x (q,u). Here, we 
set h = 0, and the system size is N = 32. 


distributed to a high energy region. Thus the distribution of the spectral weight is asymmetric 
with respect to qo mm ■ Note that the intensity extends to high energies beyond the upper 
boundary of the two-spinon continuum for q > qo, which is indicative of a four-spinon continuum 
suggested in LiCuV 04 [?]. 

In Figs. (He) and[U(d), we show S z (q,uj ) and S x (q,u), respectively, at A = 0.6 for the XY-like 
case, where the system is in a gapless chiral phase nu. Both S z (q,LU ) and S x {q,oj ) resemble 
those in the isotropic Heisenberg case in their overall spectral shapes, while they move to lower 
energy. We observe that the intensity at the lowest-energy peak is suppressed for S z (q,uj), while 
it is enhanced for S x (q,u). Note that S z (q,uj ) and S x (q,uj ) have the same total intensity, and 
the difference of the distribution of the spectral weight in the ( q , oj) space is essential. The spin 
excitation spectra are supposed to be gapless in the gapless chiral phase, but we do not clearly 
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Figure 2. The dynamical spin structure factors S z (q, oj) and S x (q,u ) at J\ = —1, J 2 = 1, 
and m = 0.25 for typical values of A: (a),(b) A = 1 for the isotropic Heisenberg case; (c), (d) 
A = 0.6 for the XY-like case; and (e), (f) A = 1.4 for the Ising-like case. Three panels in the 
left column [(a), (c), (e)] are S z (q,uj), and those in the right column [(b), (d), (f)] are S x (q,uj). 
Here, we set h = 1.023 for A = 1, h = 0.887 for A = 0.6, and h = 1.177 for A = 1.4, and the 
system size is N = 32. 


see the gapless nature in the present results because of the finite-size effect. In Figs. [D(e) and 
m, we show S z (q,u) and S x (q,u), respectively, at A = 1.4 for the Ising-like case, where the 
system is in the so-called uudd phase nans]. The overall spectral shapes are again similar to 
those of the isotropic Heisenberg case, while they are pushed up to higher energy in contrast to 
the XY-like case. In addition, the intensity at the lowest-energy peak is enhanced for S z (q,u), 
while it is suppressed for S x (q, u>), which is the reverse of what we find in the XY-like case. In 
the uudd phase, the longitudinal spin excitation spectrum is gapless, while the transverse one is 
gapped, but these features are not clear because of the finite-size effect. The detailed finite-size 
scaling analysis is an important future problem, which will be reported elsewhere. 

Now, we move on to the investigation in the magnetic field. In Fig. [21 we present intensity 

























Figure 3. The spin excitation energy Ai(N,M) at J\ = — 1 and J 2 = 1 for several values of 
A. Here, we show the results only in the field-induced SDW/SN regime. Ai (IV, M) in the VC 
phase is supposed to be zero in the thermodynamic limit. The system size is N = 128. 


plots of S z (q,u ) and S x (q,u ) at to = 0.25, corresponding to the field-induced SDW regime. 
Figure [2j a) shows S z {q,ui ) at A = 1 for the isotropic Heisenberg case. We find a lowest-energy 
peak at {q o,u>o) = (97r/33,0.00). q 0 moves toward small momentum from the position at zero 
magnetic field. By analyzing q$ for various values of m, we have confirmed that qo changes with 
m following a relation qo = (1/2 — m )it in the field-induced SDW/SN phase [8], which supports 
the bosonization result [2|. On the other hand, ujo is nearly zero, indicating a gapless mode of 
the longitudinal spin excitation. Note that we also find a ferromagnetic peak at the origin due 
to the finite magnetization. In Fig. [2](b), we show S x (q,u >) at A = 1. We find a lowest-energy 
peak at (qo,wo) = (167r/33,0.22). In particular, we have verified that uq is extrapolated to a 
finite excitation energy in the thermodynamic limit, indicating a gapped mode of the transverse 
spin excitation. We notice that S x (q,u ) = [S + (q,uo) + S~(q, w)]/4, i.e., it represents the spin 
excitation changing the magnetization by one. Since a finite energy is required to break a bound 
magnon pair, we have an excitation gap corresponding to the binding energy. 

All these features found in the isotropic Heisenberg case are similarly observed for both the 
XY-like [Figs. [3(c) and [Sfd)] and Ising-like [Figs. [2](e) and [51(f)] cases. That is, for S z (q,u), qo 
changes with m as qo = (1/2 — m)7r in the field-induced SDW/SN phase, and ujq is nearly zero, 
indicating a gapless mode. For S x (q, w), ojo is found at a finite energy, while it is extrapolated 
to a finite excitation energy in the thermodynamic limit, indicating a gapped mode. Here, we 
clearly see the difference of ujq of S x (q, cu) between the XY-like and Ising-like cases. Compared 
to the isotropic Heisenberg case [Fig. [2Kb)], ujq decreases in the XY-like case [Fig. (2)(d)], while it 
increases in the Ising-like case [Fig. 0(f)]. Thus the excitation gap is suppressed in the XY-like 
case, and it is enhanced in the Ising-like case. 

To clarify the dependence of the excitation gap of the transverse spin excitation, we calculate 
a spin excitation energy, defined by 

Ai(IV, M) = [E 0 {N,M + 1) + Eq(N, M - 1) - 2E 0 {N,M))/2, (4) 

where Eq(N, M) is the lowest energy of the IV-site system in the subspace of M at h = 0. Here, 
we use an ordinary DMRG method, not dynamical one, to obtain Ai(IV, M). In Fig. [3l we plot 
Ai (N,M) for several values of A. The results in the field-induced SDW/SN regime are shown. 
Note that we find reentrant behavior in the XY-like case, i.e., the ground state changes from the 
VC state to the SDW/SN state and turns to the VC state again as the magnetization increases 
from zero to the saturation [5]- We find that Ai(V, M) increases as A increases over the XY-like 











and Ising-like regions. This indicates that the easy-plane anisotropy disfavors the formation of 
two-magnon bound states, and the easy-axis anisotropy stabilizes the two-magnon bound state, 
which is consistent with the behavior of correlation functions [5]. 

4. Summary 

We have investigated the spin excitation dynamics of the spin-1/2 one-dimensional J\-Ji XXZ 
model in the magnetic field by numerical methods. We have observed that in the field-induced 
SDW/SN phase, the longitudinal spin excitation is gapless, while the transverse spin excitation 
is gapped, irrespective of the exchange anisotropy. The excitation gap of the transverse spin 
excitation increases as the exchange anisotropy increases, which is indicative that the easy-axis 
anisotropy stabilizes the two-magnon bound state. 
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